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Abstract. For two particles in a disordered chain of length L with on-site interaction U, a duality trans- 
formation maps the behavior at weak interaction onto the behavior at strong interaction. Around the fixed 
point of this transformation, the interaction yields a maximum mixing of the one body states. When L « Li 
(the one particle localization length), this mixing results in weak chaos accompanied by multifractal wave 
functions and critical spectral statistics, as in the one particle problem at the mobility edge or in certain 
pseudo-integrable billiards. In one dimension, a local interaction can only yield this weak chaos but can 
never drive the two particle system to full chaos with Wigner-Dyson statistics. 
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The competition between two body (electron-electron) 
interaction and one body kinetic energy in disordered sys- 
tems is a fundamental problem of permanent interest. We 
denote by U, t and W the parameters characterizing the 
interaction, the one body kinetic energy and the fluctu- 
ations of the random potential for a d-dimensional sys- 
tem of size L. When U is small, the TV-body eigenstates 
are close to the symmetrized products of one body states 
(Slater determinants for spinless fermions) which contain 
the effects of t and W completely. The effect of U can be 
treated as a perturbation, yielding a mixing of those sym- 
metrized products. When U increases, the consequence 
of this mixing is that an increasing number of one body 
states is needed to describe the exact A^-body states. If 
the one body states are localized by the disorder, delo- 
calization in real space results from this mixing. This is 
why the interaction can induce in certain cases metallic 
behavior in a system which would be an insulator other- 
wise. When U is large and dominates, one can get on the 
contrary a correlated insulator which might be metallic at 
weaker interaction. A Wigner crystal pinned by disorder is 
a good example of such an interaction-induced insulator. 
In the large-?/ limit, t becomes the small parameter, and 
one expands in powers of /U . The issue is to know the 
range of validity of these perturbative approaches, and to 
describe how the system goes from the first limit to the 
second when U increases. For this purpose, we consider 
a one-dimensional disordered lattice with on-site interac- 
tions. 
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We discuss the simple case of two electrons with oppo- 
site spins (the orbital part of the wavefunction is symmet- 
ric as in the case of two bosons). For the main sub-band 
of states centered around E — 0, both in the limits where 
[/ = (free bosons) and U — co (hard-core bosons), the 
two body states can be described in terms of two one body 
states. We use a duality transformation U at^ /U to 
map the small JJ-limit onto the large tJ-limit (a w %/24). 
We first prove that the lifetimes of the free boson states 
and of the hard-core boson states are equal at the fixed 
point Uc of the duality transformation. At Uc one has the 
maximum mixing of the one body states by the interaction 
and the enhancement factor [1] is maximum for the two 
particle localization length L2. Far from Uc, L2 is smaller 
and satisfies the duality relation. L2 — > Li (the one par- 
ticle localization length) both for [7^0 and U ^ 00. 
The study of the signature of this duality transforma- 
tion on the spectral fiuctuations is very interesting. For 
E — taking L ~ Li and increasing J7, one gets two 
thresholds defining a range of interaction Up < U < Ua. 
Outside this range, the levels are almost uncorrelated. In- 
side this range, the level repulsion is maximum, but does 
not reach the universal Wigner-Dyson (W-D) repulsion. 
The two particle system is not fully chaotic, but exhibits 
a weak chaos which is not arbitrarily situated between 
Poisson (integrable) and Wigner (chaos) . The spacing dis- 
tribution p{s) between consecutive energy levels and the 
statistics S2 [E) (variance of the number of energy levels 
inside an energy interval E) are characteristic of the third 
known universality class [2]. One finds p(s) « 4sexp(— 2s) 
and S2{E) ~ 0.16-1- 0.41i? for periodic boundary con- 
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ditions. This is very close, if not identical, to the distri- 
butions found in many "critical" one body systems, such 
as an electron in a 3d random potential at the mobility 
edge [2,3,4] or in certain pseudo-integrablc quantum bil- 
liards (rational triangles [5] , rough billiards [6] and Kepler 
problem [7]). Furthermore, p{s) saturates to 4sexp(— 2s) 
for U ~ Uc only when the ratio 1 < Li/L < 10. We 
conclude that a local interaction can never drive the two 
particle system to full quantum chaos with Wigner-Dyson 
statistics in one dimension, but can at most yield weak 
critical chaos in a certain domain of interaction and of the 
ratios Li/L. We show in addition that this weak chaos is 
accompanied by multifractal wavefunctions, in agreement 
with Ref. [8]. 

Each one particle Hamiltonian is given by 

L-l L 

Ho = tJ2i\n){n + 1| + |n + l)(n|) + ^ K|n)(n| (1) 



n=l 



where Vn is uniformly distributed inside [— Wi+M^] and 
the interaction is described by 



-ffint = XI \nn){nn\ 



(2) 



Denoting \a) the one body state of energy and the 
amplitude (n\a) = ^a{n), only two one body states \a) 
and 1/3) are necessary to describe a free boson state \fh) = 
I a/3) with components (nin2|/6) given by 



(3) 



In this free boson basis, the interaction matrix elements 



(a/3|ifi„t|7^) = 2[/Q2:l where 



(4) 



n=l 



This Q-matrix has been studied [8] in one dimension for 
L > Li. It was found that its lines, where are the terms 
*5a^ coupling a given |a/3) to all the other |7^), is a mul- 
tifractal measure in the \ j6) space. Therefore the effective 
density of \fb) states coupled by the square of the in- 
teraction matrix elements is reduced. For L « Li, one has 
p2^{Li) oc Li"^^, instead of the total two body density 
p2 oc Lf . This was confirmed by a study of their Fermi 
golden rule decay. 

As noticed in Ref. [9], there is a useful representation 
for an on-site interaction — > ±oo, composed by a small 
set of Lm = L "molecular" states \nn) and by a large 
set of Lh — L{L — l)/2 hard core boson states \hc) built 
from re-symmetrized Slater determinants. Their compo- 
nents {nin2\hc) are given by 



{^a{n2)^l3{ni) - ^a{ni)^f3{n2))sgn{n2 - m) 
V2 



(5) 



The re-symmetrization is insured by the function 
sgn(a;) := a;/|a;|. For hard wall boundaries (which we as- 
sume in our discussion), one has to use the same one body 



states \a) and |/3) both for \ fb) and \hc). However, if the 
system is closed on itself and forms a ring threaded by 
a flux the re-symmetrization on a torus reveals inter- 
esting topological aspects (see ref. [10]). One has to use 
different one body states in the two limits which differ in 
by half a flux quantum (to periodic boundary conditions 
for [/ = correspond anti-periodic boundary conditions 
for U — » ±oo). 

In this basis, the two body Hamiltonian H has the 
structure 

Hm He 



n 



Hq Hh 



(6) 



Hm and Hn are Lm x Lm and Ljj x Lh diagonal matrices 
with entries U + 214, and ehc = + e/3, respectively. Hm 
and are coupled by a rectangular matrix He result- 
ing from the kinetic terms i X^^Z^ (|n, n + l)(n,n| -|- 
|n, n){n,n+ 1| -|- h.c.) of H. The matrix elements of Hq 

are equal to \/2t (^il'i3{n)D^a{n) - ^a{n)D^/3{n)^ , where 

D'^a{n) ■■= \i'ain+l)-^a{n-l). The states \hc) of energy 
Efic ~ are coupled by a term of order t to the few states 
\nn) of energy w U. Their lifetime becomes infinite when 
U ±oo. 

Projecting an eigenstate \A) of energy Ea onto the 
states \nn) and \hc) 



\A) = J2c'X\nn)+Y,c\^\hc). 



(7) 



71=1 hc=l 

one finds the relation 

'Hm + Jm{Ea) 

Hu + MEa) 



'-'A 



Ea 



which holds for arbitrary U. and are vectors of 

Lm components and Ljj components c'^ respectively. 
The matrix Jm{Ea) has L^ elements given by 

T /T-, N v-^ ((nn\Hc\hc)(hc\Hc\mm)) 
JM{EA,nn,mm) = - — ' J ' ^ (8) 



hc=l 



Ea — e-hc 



and the matrix Jh(Sa) has L^ elements of the form 

TIT? h h\ ^ {hc\Hc\nn){nn\Hc\hc) 
MEaMM = Y U + 2Vn-EA • 

n=l 

Wc consider the main sub-band of states with energy 
Ea « 0, resulting from the mixing by a perturbation of 
order t^/U of a few states \hc) for which Jyi{Ea) can be 
simplified. Assuming W,t <^U and U +2Vn — EA « U, one 
has to evaluate J2n=ii^^\-^c\nn) {nn\Hc\hc) . This expres- 
sion is composed of 12 sums over n, each of them having 
the form 



L 

E 

71=1 



<I^„{n)^0{n')^^{n)Mn") 



(10) 



with various combinations of n',n" = n ± 1. Therefore 
(52^^ is not exactly equal to the Q^^^ occuring around the 



Xavier Waintal et al.: TIP II: Critical statistics and maximum mixing of the one body states 



3 




Fig. 1. Fq (open symbols) and .Too (filled symbols) for L = 25 
(triangles), L = 76 (diamonds), L = 150 (squares) and L = 
200 (circles). Dashed (dotted) lines are the Fermi golden rule 
expression of Eq.l3 (Eq.l4). 



free boson limit. However, this difference should not be 
statistically relevant and if one neglects it, one finds 



Jn{EA,hc, he) « ±2(24) 



1/2 



(11) 



This establishes the duality transformation U 
\/2A{f' /U) for ft: which maps the distribution of the 
coupling terms between the when U is small onto the 
distribution of the coupling terms between the \hc) when 
U is large. To illustrate this duality around the fixed point 
Uc = (24) ^/"^t, we have nmncrically calculated the average 
over the disorder of the local density of states 

Pa{E) = f HE + Ea- 6^-60) (12) 

a/3 

for L — Li and Ea ~ 0. We have observed that {pa{E)) 
can be fitted by a Lorcntzian curve of width Fq (if a/3 
denotes the states |/6)) or F^o (if ctfi denotes the states 
\hc)). Fermi's golden rule gives for the widths F a.t L Li 



Fo « 27r[/2 J^pf (Li) 



2tt- 



24r 1 



pf{Li 



(13) 
(14) 



where pf^ « Ll^'^^/t ^ P2{Li) » L\/t because of 
the multifractality of the Q-matrix. As shown in fig- 
ure 1, the widths F obey the duality property ro{U) = 
F^^-s/OAt^ /U), and are described by the above golden rule 
approximations. When U ~Uc, the lifetime of the free bo- 
son states is equal to the one of the hard core boson states. 



When U increases, the statistical ensemble associated 
to the two particle system exhibits a crossover from one 
preferential basis (the free boson basis) to another prefer- 
ential basis (the hard core boson basis) . At ?7c , the statisti- 
cal ensemble is in the middle between the two preferential 
bases, the mixing of the one body states is maximum, and 
the localization length L2 therefore is maximum. When 
E = {) and U varies, a transition occurs in the two body 
spectrum in Id which is somewhat reminiscent of the one 
body case in 3d when E = Q and W varies. In the two 
cases, one can expand in powers of U/t or t/U , of W/t or 
t/W when the system is under or above the critical thresh- 
old, respectively. The analogy is of course not complete, 
since the two body spectrum in Id has essentially Pois- 
son statistics in the two perturbative limits, while the one 
body spectrum in 3d exhibits W-D statistics iov W <^ t 
and Poisson statistics when W ^ t. Nevertheless, the 
question whether or not the spectral fluctuations are of 
the same kind in the vicinity of the threshold deserves to 
be investigated. 

Before showing the results, two arguments can be men- 
tioned: (i) A multifractal Q-matrix is incompatible with 
Wigner-Dyson level repulsion in Id. The states \fh) or \hc) 
are directly coupled by U or /U to an effective density 
< p2- Therefore, nearest neighbors in energy are very 
likely uncorrelated, enhancing the probability to find level 
spacings small compared to their average, (ii) Looking at 
Eq. (2), one may assume that a broad distribution of the 
matrix elements of Hy^ may yield a single dominant cou- 
pling term. This allows us to consider that the states \hc) 
are mainly coupled via a few states \nn). This is not far 
from the case discussed by Bohr and Mottelson [11] (cou- 
pling via a single state) where the consecutive levels Ea 
of H satisfy the conditions euc < Ea < ehc+i < Ea+i, 
Che being the consecutive levels of Hu- Since the statistics 
of Hu is essentially Poissonian, this forbids to have the 
Wigner-Dyson rigidity for the Ea- The most rigid spec- 
trum would be achieved by putting the Ea exactly in the 
middle of two consecutive Chc- It is straightforward to find 
thatp(s) would then be equal to pc(s) = 4sexp(— 2s). This 
"semi- Poisson" distribution, where Pc('S) oc s for .s <C 1 (as 
the Wigner surmise Pw{s) = (7r/2)s exp(— 7r.s^/4)), and 
which decays as exp(— s) for s » 1 (as the Poisson dis- 
tribution pp{s)), is attracting increasing interest since it 
characterizes [5] weak chaos in systems which arc "crit- 
ical" (i. e. which can be mapped [7] in a certain "inte- 
grable" basis onto an Anderson model at its critical point). 
It was shown in Ref. [12] that a plasma model for the 
energy levels can reproduce the level statistics of a dis- 
ordered conductor if the logarithmic Wigner-Dyson level 
repulsion was screened at an energy scale given by the 
Thouless energy Et- Since Et is of order of the mean 
level spacing Ai at the mobility edge, one can expect that 
a short range plasma model should also give this "semi- 
Poisson" distribution. This was recently confirmed [5] for 
a logarithmic level repulsion truncated to the first con- 
secutive level, where one finds p{s) — 4sexp(— 2s) and 
S2{E) = 1/16 + E/2- After averaging over certain bound- 
ary conditions, Pc{s) was also found to give [4] a good fit 
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Fig. 2. TIP-p(s) for L = Li = 150 and U = 1: Continu- 
ous, dotted and dashed lines correspond to the "Semi-Poisson" , 
Wigner and Poisson distributions respectively. Full (empty) 
circles correspond to hardwall (periodic) boundary conditions. 
Inset: TIP-p(s) (circles U — 1) compared to the one particle 
p{s) (squares) for L = Li = 150. 
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Fig. 3. E2{E) for L = Li = 100 and U = 1.25 E in units 
of the mean level spacing. Full (empty) circles correspond to 
hardwall (periodic) boundary conditions. 

for the p{s) of the one body spectrum in 3d at the mobility 
edge. 

These observations lead us to study p{s) as a func- 
tion of U and of the ratio Li/L around E = 0, where p2 
has a van Hove singularity ioi W = U = 0. The disorder 
and the interaction remove the singularity and p2 devel- 
ops a small plateau around « 0. It is inside this plateau 



that p{s) and 2J2 have been investigated. For U — 1 and 
L ~ Li ~ 150, the spacing distribution p{s) is shown in 
Fig. 2, in good agreement with the "semi-Poisson" distri- 
bution. In the inset, one can see that this is specific to the 
two particle problem, and does not characterize the single 
particle spectrum for the same chain of size Li. 

The statistics S2{E) displays also the behavior that 
one expects at a critical point or for a short range plasma 
model. E2{E) « + XiE is shown in fig. 3. The val- 
ues of the compressibility xi = '^^"che^oo ^2{E) / E seem 
to depend on the boundary conditions for the first con- 
secutive levels, as for the Anderson transition in 3d. The 
values for the two particle system 0.41 (periodic) and 0.45 
(hard wall) coincide with the values given [4] by Braun et 
al. for certain combinations of boundary conditions in the 
different directions. This sensitivity of xi to the bound- 
ary conditions may disappear for larger E (smaller time 
scales), as it is the case for the one particle spectrum at 
the mobility edge [13]. 

We now study the domain of validity for weak chaos 
and universal critical statistics. To measure the deviation 
of p(s) from the -Pp(s) or Pvv(s)j we use the functional 

/o ds(pp(s) -pw(s)) 

with b = 0.4729. A Poisson spectrum gives 77 = 1, a 
Wigner-Dyson spectrum gives 77 = and "semi-Poisson" 
corresponds to rjc = 0.386. 

First, we vary J7, imposing the relation L — Li. We 
assume Li to be given by the weak disorder formula 
Li « 25/W'^. In the hmits C/ ^ and C/ ^ cx), the 
TIP-levels are given by £(,/ ~ euc — + e/3 which are 
uncorrelated (at least on energy scales smaller than the 
one particle level spacing Ai). This gives pp(s) for the 
TIP-spectrum. One can see in fig. 4 that p{s) deviates for 
small U and /U from pp(s) observed at J7 = and 
U = 00. rj{U) exhibits a plateau (C/p < U < ?7h) at 
the value r^c » 0.386 which characterizes Pc{s))- We sug- 
gest that U-p and J7h are given by the conditions which 
hold [14] for systems in which there is no coupling term 
between consecutive energy levels. In agreement with the 
general picture developed in Ref. [14] (see also Ref. [15]), 
the threshold appears when the strength of the coupling 
terms becomes of the order of the spacing 1/ p'^ of the 

directly coupled levels. This gives 2Uv/lI''^ « U/L\ '^^ 
and 2^/2it^/{U^L\''^) w At/L\'^^, respectively, and seems 
to account for the size of the plateau of the curve r?(t/). 
This favours a line of critical points rather than an iso- 
lated point when the condition L = Li is imposed. The 
[/-dependence of xi exhibits the same plateau as ri{U). 
The results shown in fig. 3 for C/ = 1.25 do not vary inside 
the plateau. 

Second, for U — 1.25 and L = 100, we vary the disor- 
der parameter W in order to study the role of the ratio 
Li/L. When L ^ Li, the small fraction of TIP-states re- 
organized by U disappears behind the set of TIP-states 
which are not re-organized by U and which remain uncor- 
related. p{s) does not differ from pp{s). On the contrary. 
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Fig. 4. Weak chcios and interaction for L = Li: Crosses, pluses 
and stars represent rj{U) for L = 100, L = 150 and L = 200, 
respectively. 
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Fig. 6. d{q) at f/ = 2 for L = Li = 120 in the bases \ fb) (open 
symbols) and \hc) (filled symbols). The inset shows the estab- 
lishment of the power-law behaviour of l2{D) for increasing 
values of U in the \hc) basis. 




Fig. 5. Weak chaos and disorder for U = 1.25: r]{W) as a 
function of Li /L with L = 100. 



when i <C Li, we approach the clean hmit which can be 
easily understood: For W = 0, the total TIP-momentum 
K is conserved and the Hamiltonian matrix has a block 
diagonal form in the free boson eigenbasis. Denoting by 
e-y = 2 cos k-y the eigenenergy of a single particle with mo- 
mentum k-y, and assuming periodic boundary conditions, 
one can easily show that the {L + l)/2 TIP eigenenergies 
Ea{K) of total momentum K modulo (27r) satisfy 



2U 



E 



= 1. 



(16) 



The summation is restricted to the momenta kj + kg = K 
modulo (27r). This sequence of {L + l)/2 levels exhibits 
the same property than in the case discussed by Bohr and 
Mottelson: For a given K, the levels for U alternate 
with the levels for [7 = 0. This sequence oi Ea{K) can only 
exhibit level repulsion between next nearest neighbors, but 
not the long-range repulsion necessary for Wigner-Dyson 
statistics. Since the total spectrum is the imcorrclated su- 
perposition of the L different scqiiences associated to dif- 
ferent total momenta K, one gets for the TIP spectrum 
uncorrelatcd levels for any value of U when W = 0. 

Between the two limits L ^ Li and L ^ Li, one 
can see in fig. 4 a plateau where the spectrum is more 
rigid, rj{Li/L) saturating to the "semi-Poisson" value. 
This plateau for U = 1.25 is observed for 1 < Li/L < 10. 

A critical statistics for the spectrum is usually related 
to multifractal wave functions. We have studied the pro- 
jection C^^ of an eigenstate \A) onto the states \fb) and 
onto the states \hc). Since the perturbation matrix ele- 
ments (divided by C/ or i^/[/) are multifractal, one expects 
that this will eventually yield a multifractal structure for 
\A) in the two bases for large and small enough U, respec- 
tively. We have done the same analysis as in Ref. [8] to 
which we refer for technical details and references. We di- 
vide the plane (a, (3) into {L/D)'^ boxes of size D and com- 
pute for each box i the probability Pi = X]a/3Gboxi l^'^^l- 
After ensemble averaging, the moments Iq{D) = X^iP? 
have a power-law behaviour if U is large enough for the 
1/5) basis (inset of fig. 3), and if U is small enough for the 
\hc) basis. So one finds Ig{D) oc £)(9-i)<^(9), d{q) defin- 
ing the corresponding generalized Renyi dimensions. The 
d{q) are roughly equal in the two bases (fig. 3) and do 
not depend on U as soon as the power-law behaviour is 
established for Iq{D), i. e. inside the line of fixed points. 
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In summary, wc have shown that this simple two par- 
ticle system in Id is characterized by two perturbative 
regimes (in U or /U). These regimes are dominated by 
the existence of a single preferential basis and arc related 
to each other by a duality transformation. In the mid- 
dle, the mixing of the one body states by the interaction 
is maximum, but the system is not chaotic (no Wigner- 
Dyson statistics) . It exhibits the weak chaos of the critical 
systems, as confirmed by the spectral fluctuations. The 
curves r}{U,W),xi{U,W) and dq{U,W) in the \fh) and 
\hc) bases are consistent with the existence of a weak crit- 
ical chaos for a given size in {U, W) domain located around 
Uc and Li/L « 1. In one dimension, a local interaction 
can never drive a two particle system towards full chaos 
with Wigner-Dyson statistics. 

We acknowledge useful discussions with E. Bogomolny, G. 
Montambaux and P. Silvestrov. 
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